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We suggest that Free Random Variables, represented here by large random matri- 
ces with spectral Levy disorder, may be relevant for several problems related to 
the modeling of financial systems. In particular, we consider a financial covariance 
matrix composed of asymmetric and free random Levy matrices. We derive an al- 
■ gebraic equation for the resolvent and solve it to extract the spectral density. The 

, free eigenvalue spectrum is in remarkable agreement with the one obtained from the 

Q [ covariance matrix of the SP500 financial market. 
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5r , Noise with power law distribution is encountered in numerous stochastic sys- 

tems in physics, biophysics and finances (e.g. financial time series) [1,2]. One 
of the crucial problems in such systems is an understanding of the character of 
correlations between numerous agents. It is therefore tempting to try to find 
a description, in which we could follow the stochastic evolution not of a single 
object, but rather an evolution of an array (matrix) built of many objects. 
Such description requires basically a generalization of the classical concepts of 
theory of probability to non-commuting variables. This fundamental concept 
has been introduced by Voiculescu, using the powerful theory of free ran- 
dom variables [3]. Another major twist happened when it was realized, that 
abstract non-commuting operators representing FRV have an explicit realiza- 
tion in terms of certain large random matrices, when the size of the matrix 



^ Talk at NATO Advanced Research Workshop "Applications of Physics to Eco- 
nomic Modeling", Prague, 8-10 February, 2001. 



Preprint submitted to Elsevier Preprint 



1 February 2008 



tends to infinity. Free random matrices with Gaussian fixed point, have been 
appUed to many physical problems [4-6]. Recently, using some of the concepts 
developed in [3], we extended the concept of free random Levy variables to 
matrices [7] , and suggested that the results may be relevant for addressing the 
issue of noise in stochastic systems with power law distributions. In this talk, 
we explain in a pedagogical way what is the corner-stone of stability for non- 
commuting variables. Then we move towards a sample application, discussing 
the financial correlation matrix. 

Let us consider a Gaussian distribution in classical theory of probability, i.e. 
the distribution 



P(^) = ;y|^e ^ (1) 



with the variance a. Now let us ask the question, what is the distribution of 
the sum of two identical independent Gaussian distributions, i.e. what is the 
distribution of 



^2 = a;i + X2 (2) 

where P2('S'2) = V^^i) 'S>p{x2)- Since the Fourier transform unwinds the con- 
volution, if we introduce the characteristic functions 



P{x) = ^ / 4>{q)e-''''dq , (3) 
then, $2(9) = 0(5) ■ The inverse Fourier transform gives 



i.e. again the normal distribution but with the new variance (T2 — V^cr. One 
may say, that the logarithm of the characteristic function 



In $2 (g) = In 0(g) -Fin 0(g) (5) 

realizes a certain additivity law for Gaussian variables (2). 

Is it possible to find an analog of this construction for random matrix ensem- 
bles? Consider again a Gaussian ensemble built of large N x N matrices M, 
for which the pertinent resolvent is generically given by 
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G(.) = / dMe-"""'! 1, (-^) , (6) 

where the the potential V{M) equals ^TrM^. For simplicity we put the scale of 
the Gaussian distribution equal to one. This Green's function fulfills a simple 
algebraic equation G{z) = {z — G{z))~^, and the normalizable solution reads 
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Giz) = -iz-^/^^). (7) 

Using the fact, that the averaged spectral distribution is related to the dis- 
continuities of the Green's function^ 



p(A) = -- lim lmG{z)\z=\+ie 
we recover Wigner's semicircle law 



p(A) = -v/4^. (9) 

Let us now try to add two independent matrix Gaussian ensembles, as we 
did in the case of commuting variables of the classical probability theory. In 
other words, we are asking what is the distribution of the eigenvalues, if the 
corresponding Green's function is given by 



NVi(Mi) -NV2{M2) 



— Tr f -I . (10) 

This time the "convolution" is matrix- valued, and the matrices Mi, M2 do 
not commute! 



Free Random Variable (FRV) calculus provides a surprisingly simple solution 
to this problem. First, let us find the functional inverse of the Green's functions 
for ensembles Mi and M2, i.e. find the solution G[B{zy\ = z for every z for 
both ensembles. Then, define R{z) = B{z) — ^jz for both ensembles. The 
function R is additive 

/2i+2W =i?iW+i?2W. (11) 
^ It is obvious using special functions property -^^^ = PV^ — i'7Td{x). 
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Proceeding now in the reverse order one can reconstruct 6*1+2(2;). 



We exemplify this construction considering two identical Gaussian ensembles 
(6). Since for both the Green's function is given by (7), we get 

Bi{z) = B2{z) = z + l/z (12) 
so Ri{z) — R2{z) — z and the additivity law gives 

Ri+2{z) ^ Ri{z) + R2{z) ^ z + z ^ 2z (13) 

or, equivalently, Bi^2{z) = 2z + 1/ z. By substitution z — > Gi+2(^) and by the 
definition of the functional inverse, we immediately get 

G,+2{z)^\{z-V^^). (14) 
The discontinuity of this function yields the spectral distribution 



Pi+2(A) = -y8^. (15) 

The matrix convolution of two Gaussian ensembles is again a Gaussian en- 
semble, but rescaled by a factor \/2. 

The additivity law (11) in non-commutative probability theory of random 
matrices forms an analog of the additivity law for the logarithms of the char- 
acteristic functions in case of classical probability theory. 

This example demonstrates that the techniques of FRV offer a powerful short- 
cut when we seek the distribution of variables coming from "sum" of different 
ensembles. The above construction could be generalized to higher-order poly- 
nomial measures, sums of random and deterministic ensembles, free products 
of ensembles and also to strictly non-hermitian ensembles [5,6]. 

We focus now on the crucial question of stability. We know that in the classical 
probability theory the concept of Gaussian stability (central limit theorem) 
was generalized by Levy and Khintchine to distributions with a power-law tail. 
Is it possible that the Gaussian stability for FRV allows for a similar general- 
ization? Works of Bercovici and Voiculescu [3] provide a positive, but formal 
answer to this problem. Recently, we have suggested an explicit construction 
for the matrix measure of such ensembles, borrowing on the Coulomb gas anal- 
ogy well-known in random matrix theory. Generically, the potential is given 

by 
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V'(X) = 2 Re ^(A + iO) 



(16) 



and for free Levy matrices, G{z) satisfies an algebraic equation in large N 
limit [7] {a 1) 



b C'iz) -{z-a) G{z) + 1 = 0, (17) 

in the upper half-plane, and follows by Cauchy reflection in the lower half- 
plane. The parameter b is related to the Levy index a, asymmetry f3, and 
range or scale 7. Further details are discussed in [7]. We would like to mention, 
that the explicit form of the matrix measure is known only in few cases. 
Even in these cases. Levy ensembles exhibit very nontrivial behavior, and the 
potentials are non-analytic, contrary to the broadly studied polynomial cases. 

Let us apply some of these observations to econophysics. Recently, it was 
pointed out that financial covariances are permeated by Gaussian noise in the 
low-lying eigenvalue region with consequences on risk assessment [8] . Here we 
show, following our recent analysis [9] , that the financial covariancc spectrum 
is throughout permeated by free Levy noise. We quantify the fiuctuations in 
a stochastic system by the use of a covariance matrix C, 



C^J = 7fL;J:MuMtJ, (18) 
t=i 

where isT x N (asymmetric). The normalization in (18) follows from the 
fact that the "variance" for Levy matrices grows like T^/" where o; = 2 is the 
expected diffusive limit (Brownian). Free Levy ensembles are super- diffusive 
for < a < 1 and sub-diffusive for 1 < a < 2 with a — 1 the critical divide. 
We define the empirical matrix of relative returns 



Mti = {xu+i - xu) /Xio - (M) , (19) 

where Xit are the raw returns made of the price x^ of stock i at time t. The 
returns are normalized to the initial price Xio to insure scale invariance. They 
carry zero mean after subtracting the average relative return {Mfj)t,i- For the 
raw prices we will use the daily quotations of = 406 stocks from the SP500 
market over the period of T + 1 = 1309 days from 01.01.1991 tiU 06.03.1996 
(ignoring dividends). For these data the matrix asymmetry is m = 1308/406 ^ 
3.22. 

The spectrum of C contains important information about the character of the 
correlations. Using Voiculescu's powerful machinery of R and S transforms 
for free random variables [3] and some of the techniques developed in our 
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Fig. 1. Left: Spectral density of FRV with a = 1/2 and different asymmetry pa- 
rameters m. Right: spectral density for several indices a, at m = 3.22. 

earlier papers we find that the resolvent for the covariance built from pure 
Levy matrices M satisfies the transcendental equation (7 = 1) [9] 



G{z) 



l + wiz) 



(20) 



where w{z) satisfies the multi- valued equation 



-e'^ •w^/"-^^ (l + w)(w + m) , (21) 

for fixed asymmetry m — T/N and Levy asymmetry parameter /3 = 0. The 
distribution of eigenvalues of the covariance matrix follows from the disconti- 
nuity of G{z). The distribution is unimodal and it's support extends to infinity. 
Since the covariance is a square, the tail distribution is characterized by an 
index a/2 since the entries Mj,, are Levy distributed with index a. 

In Fig. 1 we show some typical spectra. The case a — 2 corresponds to Gaus- 
sian distributions used recently in [8]. 

In Fig. 2 we compare the analytical results following from the free Levy covari- 
ance to the ones following from the SP500 market. Fig. 2 shows the distribution 
of eigenvalues for free Levy matrices with index a/2 = 3/4 and asymmetry 
m = 2, versus the raw SP500 data (left) and the reshuffled SP500 data (right). 
The reshuffled data are obtained by randomly permuting the time ordering 
of the price series, independently for all stocks, destroying all inter-stock cor- 
relations. Our optimal fit preserves the index of the cumulative distribution, 
albeit for a smaller asymmetry (m = 2 < 3.22). It is remarkable that our free 
Levy fit to the spectrum of the covariance suggests that the time series of 
returns is power law distributed with index a = 3/2 which is consistent with 
detailed studies of the SP500 financial time series [2] . 
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Fig. 2. Raw (left) and reshuffled (right) SP500 spectrum compared to the FRV 

result with m = 2. 

In this talk, we have summarized some of our recent results on free random 
Levy matrices and finances [7,9]. A more thorough analysis and discussions 
including (non-rotationally invariant) ensembles of Levy matrices [10] can be 
found in [7,9]. We have shown that the spectral density of the SP500 financial 
covariance is in overall agreement with a free Levy distribution of index a/2 = 
3/4. Free Levy noise may be dominant in financial covarianccs, a point that 
may affect all appraisals of risk in finances. In general, we expect the results 
discussed in this talk to be of relevance to a number of problems involving 
non-commutative probabilities, whether in physics, biophysics, social sciences 
or finances. 
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